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We consider a string with fixed endpoints where the mass density and/or the elastic coefficient 
vary in a self-affine way as function of position. It is demonstrated how the eigenvalues in the 
asymptotic limit are distributed. Scaling laws for the Weyl term of the asymptotic integrated 
density of states is established and confirmed numerically. 
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I. INTRODUCTION 



An old problem in mathematical physics is how var- 
ious irregularities influence the asymptotics of the cu- 
mulative eigenvalue distribution for a physical resonator. 
This problem has many significant physical applications, 
such as wave scattering from fractal surfaces, liquid flow 
in porous media, vibrations of cracked bodies or macro- 
molecules (polymers) etc. 

In 1910 Lorentz put forward the conjecture that the 
number of eigenmodes up to some (large) frequency to, 
depends on the "volume" and not on the shape of the res- 
onator. This was later, proved by Weyl, under condition 
of a sufficiently smooth but otherwise arbitrary bound- 
aries Q] . Later Hunt et al. improved this formula by 
including a correction term which depends on lower pow- 
ers of the frequency and also on the "surface area" of the 
resonator's perimeter. 

Berry |}) , when working on wave scattering from frac- 
tal surfaces, generalized the result of Hunt et al. to fractal 
boundaries. He conjectured that for fractal boundaries 
with Hausdorff dimension Dh, the first correction term 
should scale as lo Dh . However, Lapidus has proved |J 
that the correct fractal dimension is not the Hausdorff 
dimension, but another nontrivial dimension known as 
the Minkowski dimension Dm . 

Over the last decade or so, there have been renewed 
interest in this, and related problems, both in the math- 
ematics PHlC| and physics |Tl]~[rl|] communities. In this 
paper we consider a related problem where we study a 



string which has a irregular (self-affine) mass density and 
local elastic coefficient. For a physical realization, we 
could for instance think of a long vibrating ("fussy") 
polymer. 

This paper is organized as follows. In Section |1 we 
introduce some of the general theoretical background for 
the problem. We then in Section HI, make our conjecture 
for the integrated density of states (IDOS)for our self- 
affine string based on the results of Section [r| In Section 
IV we discuss the decimation technique, the numerical 
method used to calculated the IDOS. The numerical re- 
sults are presented in Section M, and the conclusion of 
this paper is drawn in Section W\ 



II. GENERAL THEORY 

In this section we review some results which will prove 
useful in the later discussion. Let T C K" be a bounded 
open set. Consider the ID elliptic differential equation 



V 2 u{x) + a(x)u} 2 u(x) = 0, 



(2.1) 



with homogeneous Dirichlet boundary conditions. The 
"weight function" a(x), will be assumed to be a positive 
real valued function. This equation has a countable se- 
quence of positive eigenvalues (eigenfrequencies) tending 
to infinity. 

Let N(uj,T) denote the integrated density of states 
(IDOS), i.e. the number of eigenmodes with eigenfre- 
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quency below ui. Lapidus and Fleckinger B have showed 
that the asymptotic behavior of N(u>,T) is 



N(w,T) ~ W(w,T) 



B n LO n 

(2?!-)" 



(a(x)) n/2 d n x, (2.2) 



as lo — > oo. Here £>„ denotes the volume of the unit ball 
in W 1 and V is the resonator domain. This term is usually 
called the "Weyl term" after Weyl who first proved this 
asymptotic behavior for a "classical" (i.e. non-fractal) 
resonator (the Weyl conjecture) 0. 



III. CONJECTURES FOR THE SELF-AFFINE 
STRING 

In this paper we will study the situation where the 
boundary is regular while the weight-function and/or the 
elastic coefficient is irregular. We will use the follow- 
ing physical model: Consider a freely vibrating string 
with fixed endpoints at x = and x = L (Dirich- 
let boundary conditions). It has an elastic coefficient 
E(x) = Eq + E%(x) and density profile p(x) = po+pi(x). 
Here E > and p n > are constants and Ei(x) and 
Pi (x) are strictly positive, self-affine functions. The equa- 
tion describing the vibrations of this "rough string" is 

-TT 4- E ( x ) 4- u ( x ) + ^ u ( x ) = ■ (3- 1 ) 
p[x) ax ax 



If E(x) = E , then p{x) = a(x) in Eq. (|2_lJ) |15| . 

We will assume that the elastic coefficient Ei(x) and 
the density p\ (x) vary in a self-affine way 20 1 . The con- 
cept of self-affinity is a scaling property. A function de- 
fined as h — h(x) is said to self-affine if it is statistically 
invariant under the transformation 



x 

h(x) 



Xx , 

\ H h(x) , 



for all positive A £ R, or equivalently 
h(x) ~ \- H h(\x) , 



(3.2a) 
(3.2b) 



(3.2c) 



where ~ is used in order to indicate statistical equality. 
The parameter H is the Hurst exponent (or roughness 
exponent). When H > 1, h(x) is not asymptotically flat 
and the surface is not statistically invariant under trans- 
lation. When H < 0, the variance of h(x) diverge when 
the interval over which it is measures goes to zero. h(x) 
is then referred to as a fractional noise. We will in our 
analysis assume < H < 1. Self- affinity is in practice 
only found within a restricted range of scales. In this 
work we explicitly introduce a lower cut-off I. The upper 
cut-off is the system size L, which is the length of the 
string. 

Asymptotically, on large scale, Ei(x) and pi{x) will 
dominate the behavior of E(x) and p{x) no matter what 
E and po are. However, we will investigate the system 



at intermediate scales where the constant terms may or 
may not dominate over the self-affine terms in E(x) and 
p(x). 

We will distinguish four cases: 

1. E Q > (Ei(x)) and po < (pi(x)}, 

2. E Q < {Ex(x)) and p < (px(x)}, 

3. E < (E-l(x)) and p > (pi(x)}, 

4. E > (Ei(x)) and p > (pi(x)}, 

where (•) denotes the averaging operator. 

We note that the left hand side of Eq. ([0]), 
(l/p)[(dE 1 /dx)(du/dx) + E(d 2 u/dx 2 )}, involves the 
derivative of a self-affine function, dE\j dx. This quantity 
scales as dE\jdx — > \ H ~ 1 dE\/ dx when x — > Ax, while 
E\ scales as E\ — » X H E\. Thus, we make the assump- 
tion that the term containing dE\/dx may be neglected 
in comparison to the term containing Ei in Eq. (3.1). 

With this assumption, we expect the IDOS for the self- 
affine string to behave as [n = 1, B\ = 2) 



N(uj,L) ~ W(lo,L) 




(3.3) 



in the asymptotic limit. Note that this expression can 
be written W(u,L) = loL/k{^/ p{x)/E(x)). Thus, the 
asymptotic behavior of the IDOS for a rough string is 
expected to behave as a classical (non-rough) string with 
a constant inverse velocity {y/p(x)/E(x)). 

We now discuss the four cases in turn. Assume now 
that we change the system size according to L — > XL. 
Then, for case (1) (Eq > (Ex(x)) and po < jpi(x))) af- 
ter a change of variable, it follows from Eqs. ( |3.2| ) and 
(O) that 



W(uj,XL) = 




d(Xx') 



X L+ ~ W(u),L) 



Using Eq. (3.3), this relation may be rewritten 
WiX- 1 -^ lu,XL) ~ W(uj,L). 



(3.4a) 



(3.4b) 



Note that Eqs. (3.4a) and ( 3.4b| ) are in principle equiv- 
alent, but open up two different possibilities for physical 
interpretation. In the former case, it is the number of 
eigenmodes which is scaled, while in the latter it is the 
frequency. 

Through similar arguments, we find for case (2) (Eq> <C 
(Ei(x)) and p < (pi (•?))), 



W(u>, XL) ~ XW(ui,L) 



(3.5a) 
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or equivalently 

W(A -1 w, XL) ~ W{w, L) . (3.5b) 
For case (3) (E <C {E 1 (x)) andpo 3> we expect 



W(w,XL) ~ X 1 ~ H ^ 2 W(u, L) , 
W(X~ 1+h/2 uj,XL) ~ . 



(3.6a) 
(3.6b) 



Case (4) (£ > (^(x)) and p > (oi (a?))) leads to the 
same behavior as Case (2), Eqs. (|3.5|). 



IV. OUTLINE OF THE NUMERICAL METHOD 



In order to do numerical simulation of our self-affinc 



string we discretize Eq. (3.1) on a lattice of N sites 



i,j = l,2,.. 



,N. 



(4.1) 



Here is the N x N matrix representation for the 
operator (1/ ' p{x)){d/dx)E{d/dx). It is tridiagonal. 

The method used in this paper to calculate the IDOS 
for the above equation is the decimation technique of 
Lambert and Weaire p6l. This method is based on a 
renormalization philosophy, where successive degrees of 
freedom are eliminated from the system, see Refs. [ ^6|Jl7| 
for more details. After removing the sites corresponding 
to k = N — M, . . . , TV the system becomes [n8[ 



A\f ) U j {u)=u 2 U i {u), 



with 



1. 



,N - M (4.2) 



AM) _ AM-X) 



AM-X) AM-X) 
■"■j^+X-M^N+X-M.j 



where the denominator is given by 



n (M) _ AM) 

u — ^iV+l-M.iV+l- 



M 



(4.3) 



(4.4) 



At the very heart of the method lies the fact that renor- 
malized system is equivalent to the original one in the 
sense that the two systems have the same eigenvalues. 
By repeating this procedure all degrees of freedom can be 
eliminated. The number of cigenmodes less then to, i.e. 
the IDOS, can now be shown to be equal to the number 
of negative denominators in the sequence {D^}fL 1 Jl9| . 
We would like to point out that this algorithm is very 
efficient. 



V. SIMULATION AND RESULTS 



In order to test the conjecture (3.4), and thus also Eqs. 
(3.5) to (|3.6|), we have calculated the IDOS for various 
system sizes, L, fixing the Hurst exponent H to the value 



0.7. We work in length units where I = 1/2 11 . Thus, the 
system size N = 2 11 = 2048 corresponds to a string of 
unit length (L = 1). We set E = 1 and (-Ei(x)) = 0.01 
for the cases where Eq 3> (Ei(x)), and Eq = 0.01 and 
(Ei(x)) = 1 for the cases where Eq <C (Ei(x)). Likewise, 
we set po = 1 and (pi{x)) = 0.01 for the cases where 
Po 3> (px(x)}, and p = 0.01 and (pi{x)) = 1 for the 
cases where po <C (pi(x)}. We averaged over 500 sam- 
ples in case (1), while we used 50 samples in the three 
other cases. 

As can be seen from Figs. 0a the IDOS is a linear 
function of the frequency to, and W(u,L) increases with 
system size L. This is consistent with the predictions of 
Eq. ( 3.4a| ). A similar linear behavior has been found for 
case (2) to (4), but these results are not shown explicitly. 

We test in Fig. [j]b the case (1) scalin g rel ation ( ft.4a|) , 
in Fig. H the case (2) scalin g rel ation ( 3.5a ), in Fig. [|a 
the case (3) scaling relati on ( |3.6a ), and in Fig. |]the case 
(4) scaling relation ( 3.5a ). 

In all cases except case (3), our numerical data are 
consistent with the theoretical predictions. In case (3), 
where Eq <C {E\{x)) and po 3> (pi(x)), data collapse was 
obtained by scaling W(ui, i)/exp(0.46L) (see Fig. ||b). 
This is very different from the predicted behavior, Eq. 
(| 3.6a ). We do not know the reason for this discrepancy 



between theory and numerical result. 



VI. SUMMARY AND CONCLUSION 



We have investigated the IDOS for a self-affine string, 
i.e. a string with self-affine variations in mass density and 
elastic coefficient. There are four relevant cases to study 
depending on whether the self-affine variations dominate 
or not in the mass density and elastic coefficient. We have 
compared numerical studies with a conjecture based on 
assuming that a result of Lapidus and Fleckinger ji] is 
valid for the self-affine string. 

We find that our conjecture works for three out of 
the four cases. However, it fails for the case where 
E < (Ei(x)) and p a » (pi(x)). 
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FIG. 1. (a) The asymptotic behavior of the integrated density of states (IDOS) when E 3> (E 1 (x)) and p <C (pi(x)) (Case 
1). (b) The scaled IDOS W(cu, L)/L 1+H/2 as function of cu. 
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FIG. 2. The scaled IDOS W(u,L)/L as function of cu for Case (2): E < (Ei(x)) and p < (pi(x)>. 
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FIG. 3. The scaled IDOS (a) W(uj, Vj/L 1 ^ 11 ^' 2 and (b) W(u>, L)/e 0A6L as function of w for Case (3): E < (E^x)) and 
Po > (pi(x)). 
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FIG. 4. The scaled IDOS W(u,L)/L as function of w for Case (4): E > (-Ei(x)) and p > (pi{x)). 
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